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Abstract. We present several models of the magnetic structure of solar coronal mass 
ejections (CMEs). First, we model CMEs as expanding force-free magnetic structures. 
While keeping the internal magnetic field structure of the stationary solutions, expansion 
leads to complicated internal velocities and rotation, while the field structures remain 
force- free. 

Second, expansion of a CME can drive resistive dissipation within the CME changing 
the ionization states of different ions. We fit in situ measurements of ion charge states to 
the resistive spheromak solutions. 

Finally, we consider magnetic field structures of fully confined stable magnetic clouds 
containing both toroidal and poloidal magnetic fields and having no surface current 
sheets. Expansion of such clouds may lead to sudden onset of reconnection events. 



1. Force-free expansing CMEs 



Expansion of magnetic clouds is one of the basic problems in space physics 
and geophysics, related to propagation of solar disturbances, the coronal 



mass ejections (CMEs) through interplanetary medium (Forbes et al. 2006 ). 



As the space craft passes through the cloud, the magnetic field strength 
is higher than average, the density is lower, the magnetic pressure inside 
the cloud greatly exceeds the ion thermal pressure and the magnetic field 
direction changes through the cloud (e.g. 



Burlaga et al. 1982). 



Spheromaks are stationary force-free configurations of plasma satisfy- 
ing the condition J = aB with spatially constant a. They are solutions of 
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Grad-Shafranov equation in spherical coordinates when the poloidal cur- 



rent is a linear function of the magnetic flux. We have found Lyutikov and 



Gourgouliatos (2011) self-similar expanding structures of axially symmet- 
ric force-free structures with spatially constant a-parameter. (This neglects 
the possible effects of plasma inertia.) For a given radial velocity of the 
spheromak's surface, parametrized by a(t), the following electric field 



a 

E = r—e r x B, 
a 



(1) 



and plasma drift velocity 



B(B e r ) -e r B 2 
v = |p r "t hi a. 



(2) 



keep the internal structure force-free. 

For the basic spheromak solutions the electromagnetic fields become 



B r = 2Boa^- cos#, Bg = 

Ba = B jx ( — ) sin (9, 
\ctoJ 

a a 
Eg = -Bs-r, Es = Bg-r, 
a a 



-B a 



Ji + arji 



sm ( 



(3) 



We stress that the fields remain force-free, J x B = 0. 

Expanding force-free spheromaks have complicated internal velocity struc- 
ture. At each point there is a well defined electromagnetic velocity v = 
E X B/£? 2 , normalized to c, Equation The solution can be parametrized 
by the velocity of the expansion of the boundary, which is v g = v a = and 
v r = R(t)dt ln(l/a) = vq, see Figures [T] Since ^ 0, the expansion of 
a spheromak induces rotation, with the different hemispheres rotating in 
the opposite direction (Figure [2]) . The flow lines are plotted in Figure [3] 
Close to the symmetry axis, velocity at small radii is directed away from 
the axis.. Total velocity at small r <C R IS 1) ~ vqt sin 9/R. The flow lines are 
orthogonal to the surface, matching the overall expansion of the spheromak. 



In the case where magnetic clouds are connected to the Sun (Chen and 



Garren 1993 ) , their topology shall resemble more an arcade, and their math- 



ematical description is more accurately approximated by a force-free mag- 



netic field confined in a torus. A number of studies (Burlaga et al. 1981 



2 



Cent. Eur. Astrophys. Bull, vol (2012) 1, 2 



MAGNETIC STRUCTURE OF CORONAL MASS EJECTIONS 




Figure T. Velocity structure of expanding spheromak: Radial velocity v r , v$, Vj, as func- 
tions of r for different 0. In the equatorial plane radial velocity increases linearly with 
radius, v r (6 = n/2) = (r/R)vo, while close to the axis velocity remains small in the bulk, 
sharply increasing to vo near the surface. Flows are antisymmetric with respect to the 
equatorial plane, vg becomes zero at r = 0.61-R 




r/R 6 

Figure 2: Angular velocity of rotation Q = (r sin 9) as function of radius and polar 
angle. 



Burlaga 1988 Lepping et al. 1990) have considered solutions in cylindrical 



geometry as an approximation to the toroidal structures of magnetic cloud. 

Magnetic clouds expand as they propagate away from the Sun. The 
simultaneous expansion and propagation lead to an increase of both radii 
of the torus. The expansions of the small and large radii need to scale with 
the same factor to conserve the magnetic flux that crosses the equator of 
the torus and the magnetic flux that crosses a section of the torus: $i oc 
B(f, : TfoRo anci ^2 of -Bz,T r 0) so if Rq is multiplied with some factor so does 

The magnetic flux conservation of the expanding structure requires that 



0, 



Cent. Eur. Astrophys. Bull, vol (2012) 1, 3 



3 



MAGNETIC STRUCTURE OF CMES 




4 



Cent. Eur. Astrophys. Bull, vol (2012) 1, 4 



MAGNETIC STRUCTURE OF CORONAL MASS EJECTIONS 



E, = Bi 



aa 



ZD . 



(5) 



2. Observational Signature 

The magnetic field measured by a magnetometer as a function of time de- 
pends on the relative velocity of the detector and spheromak boundary. 
Generally, the velocity of the detector with respect to the spheromak's cen- 
ter is dominated by the advection velocity of the spheromak with the Solar 
wind, which is typically a factor of several larger than velocity of expansion 
of the spheromak in the winds frame. 




(2) 



(1) 




Figure 4'- Left Panel: Sections of a spheromak with a meridional plane. The trajectories 
of the detector crossing the spheromak are the dashed lines. Right Panel: Magnetic field 
measured by a detector flying through an expanding spheromak. The detector enters the 
cloud at t = 0, at t = 1 it crosses the axis: curves (1), (2); or the equator: curves (3), 
(4). We investigate four cases corresponding to the numbers of the trajectories shown 
in left panel. The solid curve (1) corresponds to a detector flying on the equator of the 
spheromak and passing from the centre at t = 1. The dashed curve (2) corresponds to a 
detector that flies parallel to the equator and crosses the axis at t = 1. The dashed-dotted 
(3) curve corresponds to a detector that flies parallel to the axis and crosses the equator 
at t = 1. Finally the dotted curve (4) corresponds to a detector flying along the axis and 
passing from the centre at t — 1. In all cases the maximum magnetic field is measured at 
t < 1 before the detector crosses the middle of the cloud. 



Observations of CMEs often indicate that the maximum of B-field is 



reached before the spacecraft reaches the geometrical center (Burlaga 1995 
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Farrugia et al. 1993). It was long suspected that is a consequence of the 



expansion of the magnetic cloud while it moves past the spacecraft (Burlaga 



1995 Section 6.5.2): when the detector reaches the middle of the spheromak, 
the overall normalization of the magnetic field decrease due to expansion. 
Our model provides a quantitative description of this effect, see Figures [4] 

3. Resistive expansion 

The eruption mechanism of coronal mass ejections (CMEs) is currently an 
extremely active area of research. A variety of mechanisms have been dis- 
cussed in the literature, most of which require some form of magnetic re- 
connection during and after the eruption. The associated thermal heating 
is an important but hitherto relatively unexploited diagnostic of the erup- 
tive process, dynamics and geometry. Such heating gives rise to UV- X-ray 
emissions during the eruption, that can be remotely sensed, and also leaves 
its imprint in the charge states of ions detected in situ near 1 AU, as their 
ionization state responds to the newly heated plasma. 

In situ measurements of ion charge states can provide a unique insight 
into the heating and evolution of coronal mass ejection when tested against 



realistic non-equilibrium ionization modeling. We investigated Rakowski 



et al. (2011) the representation of the CME magnetic field as an expand- 
ing spheromak configuration, where the plasma heating is prescribed by the 
choice of anomalous resistivity and the spheromak dynamics. 

Physically, the parallel component of the electric field should be related 
to the current density through Ohm's law. Formally, the procedure described 
below breaks down the assumption of self-similarity, since the value of the 
parallel electric field component is not linearly proportional to the current 
density. Still, we assume that resistivity plays a subdominant role, so that 
the expansion remains approximately self-similar. We are interested not in 
the detailed properties of local resistive heating, but in general scaling rela- 
tions. The requirement that an expanding CME dissipates a large fraction of 
its initial energy by the time it reaches Earth orbit can be used to estimate 
anomalous resistivity. 

We chose as a test case, the 19 May 2007 CME observed by STEREO and 



ACE Rakowski et al. (2011). This model can provide the heating required 
between I.IRq and earth orbit to produce charge states observed in the 
CME flux rope. 



(i 
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Figure 5: Left: Three-dimensional rendering of the equilibrium field. Right: cross-section 
of the high-/3 equilibrium. Black lines represent poloidal field; on the right-hand-side the 
blue shading represents the toroidal field (multiplied by the cylindrical radius r sin 8) and 
on the left-hand-side the red shading represents gas pressure. 

4. Fully confined magnetic clouds 

Let us assume that a magnetic structure containing plasma and magnetic 
field in equilibrium is confined by some external pressure. This pressure 
is taken to be constant on the boundary of the cavity. Stability requires 
that they contain both toroidal and poloidal magnetic fields, while realistic 
configurations should have vanishing magnetic field on the boundary. For 
axisymmetric configurations embedded in unmagnetized plasma, the con- 
tinuity of poloidal and toroidal magnetic field components on the surface 
of the bubble then requires solving the elliptical Grad-Shafranov equation 
with both Dirichlet and Neumann boundary conditions. This leads to a 
double eigenvalue problem, relating the pressure gradients and the toroidal 
magnetic field to the radius of the bubble. We have found fully analytical 
stable solutions. This result is confirmed by numerical simulation. 

As the bubble propagates in the wind, the external pressure decreases. 
How will the internal structure of the bubble adjust to the changing external 
conditions? We find that the evolution of the cavity depends on the adiabatic 
index of the material: for T = 4/3 the cavity expands self-similar ly keeping 



Cent. Eur. Astrophys. Bull, vol (2012) 1, 7 



7 



MAGNETIC STRUCTURE OF CMES 



the ratio of the kinetic plasma pressure over magnetic pressure constant, 
while for T > 4/3 the structure will be dominated by the magnetic field 
pressure and by kinetic plasma pressure for T < 4/3. As a result, for T > 4/3 
the expansion of the cavity will lead to the spontaneous formation of electric 
current sheets, which are subject to resistive decay of the magnetic field, 
e.g. by forming large-scale reconnection layers. 



5. Discussion 

We discuss a number of models of the magnetic structure of CMEs. We 
found the self-similar solutions for expanding force-free magnetic struc- 
tures in spherical and cylindrical geometries. Under the assumption of non- 
relativistic MHD, when the dynamical effects of the resulting induced elec- 
tric field can be neglected, the structure remains force-free during expansion, 
J = aB. In case of relativistic expansion, when electric field has impor- 



tant dynamical effects (Prendergast 2005 Gourgouliatos and Lynden-Bell 



2008 ) , spheromak-type self-similar solutions may be found for constant ex- 
pansion velocity (R. D. Blandford, private communication). In contrast, 
our solutions depend on R and R measured at a given time, and thus are 
applicable to any R{t). 

In case of spheromaks, the internal velocity develops azimuthal compo- 
nent, va, 0, even though the surface expands purely radially. The fact that 
the magnetic field structure remains force-free during expansion is qualita- 



tively different from the solutions of Low (1982); Farrugia et al. (1995), 
which were limited to pure radial motion. As a result, even if at the ini- 
tial moment magnetic field is force-free, during expansion non-zero pressure 
forces appeared. 

We have found analytical solutions for magnetic cavities without surface 
currents. The cavities contain a magnetic field with poloidal and toroidal 
components and a hot plasma. The structure of the fields is such so that 
they drop gradually to zero at the end of the cavity. Because of that, there 
are no surface currents, unlike the force-free fields which require surface 
currents. The non-force- free equilibrium found here must of course have a 
plasma-/3 higher than about unity; a low-/3 equilibrium must be force-free 
in the interior with a current sheet at the boundary. 

A numerical test has confirmed the stability of the simplest (fundamen- 
tal radial mode) equilibrium. Furthermore, simulations have been used to 
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examine the behavior of the equilibrium when the gas pressure drops by a 
large factor (via cooling) and takes the bubble from the /3 <z 1 regime into 
the f3 ■< 1 regime. It is found that, as expected, the bubble becomes ap- 
proximately force-free in the bulk with a thin layer of high current density 
at the boundary and along the axis of symmetry. 
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